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ABSTRACT 


In  the  following  report,  the  authors  highlight  briefly  the  major  qualitative  and  quantitative 
features  of  the  NSWC  ocean  tide  models,  which  now  include  the  nine  components  (M2,  S2 ,  N2 . 
K2 ),  (K, ,  O, ,  P, ,  Qi ),  and  Mf  with  models  for  Mm  and  Ssa  nearing  completion.  Special  charac 
teristics  of  those  models  are  pointed  out,  which  require  indispensable  attention  in  various  appli¬ 
cations  particularly  in  coastal  areas  with  rapidly  varying  tides.  Since  the  tide  models  interpolate 
hydrodynamically  empirical  tide  data  at  over  2000  continental  and  island  stations,  certain  simple 
corrections  in  the  numerical  procedures  are  presented.  Plans  for  more  extensive  local  and  perhaps 
global  refinements  and  improvements  of  the  tide  models  are  mentioned. 

In  the  final  part  of  this  report,  the  authors  describe  the  NSWC  Global  Ocean  Tide  Data 
(GOTD)  tape  including  the  data  arrangement  and  the  standard  format.  The  report  ends  with  a 
program  that  computes  from  the  taped  harmonic  tidal  constants  instantaneous  oceanic  and 
geocentric  tides  at  specified  random  points  and  instances  in  the  Random— Point  Tide  (RPTIDE) 
program. 
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1.  INTRODUCTION 


Although  the  mathematical  modeling  of  ocean  tides  enjoyed  since  Newton  a  tricentennial 
fascinating  history  (Schwiderski  1980a),  due  to  the  measured  digital  knowledge  of  the  bathymetry 
of  the  real  ocean  basins,  the  solution  of  the  problem  was  destined  to  await  the  advent  of  modern 
digital  computers.  Fortunately,  along  with  the  necessary  computational  tool,  significantly  im¬ 
proved  numerical  techniques  and  hydrodynamical  notions  of  turbulent  dissipation  processes 
became  available  for  the  mathematically  discrete  and  realistic  definition  of  the  ocean  model 
(Schwiderski  1980a).  Nevertheless,  the  strictly  mathematical  formulation  of  the  problem  failed  to 
describe  the  true  tides  over  some  important  ocean  areas  where  low  tides  were  predicted  for  ob¬ 
served  high  tides.  The  computed  tidal  charts  conspicuously  displayed  none  of  the  long-known 
strong  distortions  and  retardations  of  ocean  tides  by  shallow  continental  shelves  (see,  e.g.,  Figure  2) 
and  narrow  ocean  ridges. 

To  overcome  those  fundamental  shortcomings,  the  author  (Schwiderski  1978a,  b,  1979a -d, 
and  1980a)  introduced  a  unique  hydrodynamical  interpolation  technique,  which  incorporates 
into  the  model,  empirical  tide  data  collected  from  harmonically  analyzed  tide  gauge  measurements 
along  shorelines  and  at  some  other  ocean  bottom  irregularities.  The  combined  mathematical- 
experimental  method  removed  essentially  all  purely  mathematical  deficiencies  and  made  the 
desired  accurate  charting  of  ocean  tides  possible  (Schwiderski  1979e,  1980b,  1981a  -  h,  j.  and 
1982a -d). 


All  earlier  investigations  started  from  the  basic  Laplace  tidal  equations  (LTEs),  which  in¬ 
cluded  only  the  astronomical  equilibrium  tide  as  its  sole  tide-forcing  potential.  Following  sug¬ 
gestions  of  Proudman  (1928)  that  there  might  exist  significant  interactions  between  the  oceanic 
and  terrestrial  tides,  Grace  (1930)  augmented  the  LTEs  by  including  simple  approximations  of  the 
terrestrial  tide,  the  terrestrial  response  to  the  oceanic  tidal  load,  and  the  corresponding  three 
gravity  perturbations  of  the  equilibrium  tide. 

Remarkably,  all  numerical  tidalists  disregarded  those  investigations  till  Farrell  (1972)  derived 
a  more  accurate  approximation  of  the  same  interaction  terms,  which  were  then  considered  by 
Hendershott  (1972),  Zahel  (1977),  Estes  (1979),  Parke  and  Hendershott  (1979)  and  others.  In 
order  to  avoid  the  involved  computation  of  the  loading  effects,  Accad  and  Pekeris  (1978)  assumed 
a  uniform  proportionality  between  the  oceanic  tide  and  its  loading  effects.  Though  this  type  of 
simplification  was  earlier  suggested  by  Takahasi  (1929)  and  Grace  (1930),  Accad  and  Pekeris 
justified  this  approach  by  displaying  the  striking  resemblance  of  the  corresponding  global  maps. 
The  author  (Schwiderski  1978a,  1979c,  d,  and  1980a)  adopted  the  same  approximation  for  all 
NSWC  ocean  tide  models,  and  Goad  (1980b)  computed  for  the  M2  ocean  tide  the  corresponding 
terrestrial  dip,  which  also  resembles  closely  its  cause. 

Since  the  first  publication  of  the  author’s  NSWC  ocean  tide  model  about  three  years  ago, 
an  unexpected  wide  range  of  applications  to  various  problems  in  oceanography,  geophysics, 
astronomy,  and  space  technology  has  been  conducted.  Fortunately,  most  of  those  applications 
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were  very  gratifying  as  they  confirmed  pointwise  and/or  globally  the  high  accuracy  and  usefulness 
of  the  computed  partial  tides  (see,  e.g.,  Goad  1980a,  b,  Melchior  et  al.  1980,  and  Wahr  and  Sasao 
1980).  Since  numerical  global  tide  models  are  naturally  limited  in  detailed  accuracy  in  certain  local 
ocean  areas,  some  foreseeeable  difficulties  emerged  particularly  in  applications  to  compute  loading 
effects  of  ocean  tides  on  the  earth  gravity  field. 

In  order  to  avoid  unnecessary  complications  in  such  applications,  relevant  qualifying  and 
limiting  features  of  the  NSWC  ocean  tide  models  will  be  pointed  out  in  Section  2.  Simple  immedi¬ 
ate  and  future  improvements  of  the  models  will  be  mentioned  and  discussed  in  Section  3.  The 
arrangement  of  the  gridded  ocean  tide  data  on  the  NSWC  magnetic  tape  is  described  in  Section  4. 
Finally  Section  5  describes  a  program  to  compute  from  the  taped  harmonic  ocean  tide  constants 
intantaneous  ocean  and  geocentric  tidal  heights  at  a  given  set  of  random  geographical  points  and 
specified  instances.  The  corresponding  user’s  guide  and  program  listing  are  presented  in  Appen¬ 
dixes  A  and  B. 
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2.  MAJOR  FEATURES  OF  T.JE  NSWC  OCEAN  TIDE  MODELS 


In  applications,  the  following  special  characteristics  of  the  NSWC  ocean  tide  models  require 
scrupulous  attention  and  evaluation.  All  these  features  are  directly  discernible  from  the  author's 
various  publications  listed  in  the  References. 

(a)  Modes  of  the  NSWC  Ocean  Tide  Models:  The  NSWC  ocean  tide  model  includes  now 
(see  Table  1 )  nine  completed  harmonic  partial  tides  of  the  semidiurnal  (M, .  S2 .  N: .  k:  i.  diurnal 
(K, ,  O, ,  P,.  0 1 )  and  long-period  (MO  species  with  models  of  the  additional  long-period  com 
ponents  (Mm,  Ssa)  nearing  completion.  All  partial  tides  are  represented  by  their  harmonic  con¬ 
stants,  that  is,  amplitudes  and  phases  which  are  recorded  gridwise  on  magnetic  tape  (Schwiderski 
1980d,  and  1981j).  Also,  gridwise  tabulated  data  charts  and  corange  and  cotidal  maps  of  all 
eleven  constituents  are  presented  in  NSWC  technical  reports  (Schwiderski  1979e.  1981a  -  g.  and 
1982a  -c). 

(b)  The  Gridded  Ocean  Basins:  For  the  recording  of  the  harmonic  tidal  constants,  the 
simply  connected  global  ocean  basins  are  divided  by  integral-degree  longitudes  and  colatitudes 
into  uniform  >°  x  1°  grid  cells.  Naturally,  the  model  excludes  all  inland  waters  such  as  the  Great 
Lakes  and  the  Black  and  Dead  Seas,  as  well  as  the  following  gridwise  disconnected  bordering 
waters:  the  Baltic,  Kattegat,  Irish,  Mediterranean,  Red.  Ceram.  Sulu.  and  Japan  Seas;  the  Gulfs 
of  California,  Persia,  Chihli,  and  lluraki;  and  the  Hudson  and  Korean  Bays.  Similarly,  many 
shallow  (less  than  5-m  depth!  and/or  narrow  waters  such  as  the  entire  Barrier  Reef  area,  the 
Gironde  Estuary,  and  the  Fjords  of  Norway  could  not  be  modeled. 

(c)  Computed  and  Empirical  Tide  Data:  The  recorded  harmonic  constants  include  all 
computed  and  hydrodynamically  interpolated  empirical  tide  data.  The  latter  values  are  specially 
labeled  in  the  printed  charts  (Schwiderski  I979e,  1981a  -  g  and  1982a  -  e)  for  accuracy  verifi¬ 
cation  and  for  proper  applications  (see  Section  3).  All  data  are  tabulated  as  representative  for  the 
centers  of  the  corresponding  grid  cells;  which  is  strictly  realistic  only  for  all  open  ocean  areas 
where  the  tides  vary  negligibly  over  the  mesh  area.  In  fact,  the  computed  data  are  rigorously 
defined  as  cell  averages,  which  may  deviate  somewhat  from  the  true  point  values  where  the  titles 
suffer  strong  distortions  and/or  retardations  over  short  distances. 

The  empirical  data  are  by  definition  (Schwiderski  1978a,  1979d,  and  1980a)  aligned 
boundary  values  at  continental  or  island  shore  stations  or  point  values  at  some  bottom  irregu¬ 
larities.  In  addition  to  their  obvious  displacement  (see  Figure  1  below),  many  of  those  data  must 
be  viewed  as  crude  averages  of  drastically  varying  high  tides.  Tliis  applies,  for  instance,  to  the 
Skagerak;  the  English  Channel;  the  Bays  of  Fundy,  Baffin,  and  Bristol;  the  Gulfs  of  St.  Lawrence, 
Aden,  Oman,  Thailand,  Tonkin,  and  Carpentaria;  the  Straits  of  Florida,  Hudson,  Formosa,  Bass, 
and  Cook;  as  well  as  most  Seas  of  the  Indonesian  Archipelago,  where  mostly  empirical  input 
data  are  listed. 
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Table  1.  Constants  of  ivtajor  Tidal  Modes 


Tidal  Mode 

A  (m) 

a  ( 1 0~  4 /sec) 

X (deg) 

Semidiurnal  Species 

M2  =  Principal  Lunar 

0.242  334 

1.405  19 

2h0  -  2s0 

S2  =  Principal  Solar 

0.1 12  841 

1 .454  44 

0 

N2  =  Elliptical  Lunar 

0.046  398 

1.378  80 

2h0  -  3s0  +  p0 

K2  =  Declination  Luni-Solar 

0.030  704 

1.458  42 

2h0 

Dirnal  Species 

K[  =  Declination  Luni-Solar 

0.141  565 

0.729  21 

h0  +90 

O,  =  Principal  Lunar 

0.100  574 

0.675  98 

h0  -  2s0  -  90 

P,  =  Principal  Solar 

0.046  843 

0.725  23 

-h0  -  90 

Q,  =  Elliptical  Lunar 

0.019  256 

0.649  59 

h0  -  3s0  -  90 

Long-Period  Species 

Mf  =  Fortnightly  Lunar 

0.041  742 

0.053  234 

2s0 

Mm  =  Monthly  Lunar 

0.022  026 

0.026  392 

s0  ~  Po 

Ssa  =  Semiannual  Solar 

0.019  446 

0.003  8921 

2h0 

A  =  amplitude  of  the  partial  tide 
o  =  frequency  of  the  partial  tide 
X  =  astronomical  argument  of  the  partial  tide. 

(h0 ,  s0 ,  p0 )  =  mean  longitudes  of  sun,  moon,  and  lunar  perigee  at  Greenwich  midnight 
h0  =  279.696  68  +  36  000.  768  930  485T  +  3.03  ■  10'4T2 
s0  =  270.434  358  +481  267.883  141  37T  -  0.001  133T2  +  1.9  •  10'6T3 
p0  =  334.329  653  +  4  069.034  032  957  5T  -  0.010  325T2  -  1.2  •  10' 5  T3 

where 

T  =  [27  392.500  528  +  1. 000  000  035  6DJ/36  525 
D  =  d  +  365(y  -  1975>  +  Int  [(y  -  1973)/4[ 
d  =  day  number  of  year  id  =  1  for  January  1) 
y  3?  1975  =  year  number 


and 

Int  ( x  |  =  integral  part  of  x 
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Ui)  Accuracy  in  Open  Oceans:  In  the  publications  (Schwiderski  (1478a.  b,  DTOa  -  e. 
1480a.  1481a  -  h.  aiul  1482a  -  il)  a  tide  prediction  accuracy  of  better  than  5  cm  anywhere  in  the 
open  oceans  has  been  claimed  for  the  dominant  M;  tide.  This  accepted  estimate  is  based  on  the 
worldwide  agreement  of  the  computed  tide  data  witli  empirical  data  which  carry  the  same 
accuracy.  Indeed,  due  to  the  somewhat  controversial  harmonic  analysis  of  tide  gauge  measure¬ 
ments.  no  higher  accuracy  can  be  claimed  even  by  the  most  recent  deep-sea  tide  data  (see.  e  g  . 
Schwiderski  1481b).  All  lesser  computed  partial  tides  (see  fable  1)  carry  relatively  the  same 
component  accuracy  as  M:.  finis,  in  superposition  a  total  tide  prediction  accuracy  ol  bettei  than 
10  cm  can  now  be  assumed  over  all  open  oceans. 

(e)  Accuracy  in  Coastal  Waters:  Due  to  the  a  priori  resolution  limits  of  the  1°  \  1  grid 
system  mentioned  under  (b)  and  (v)  above,  the  NSWC  tide  models  may  lose  some  of  their  open 
ocean  accuracy  (d)  in  coastal  waters  of  rapidly  varying  tides.  Of  course,  a  similar  loss  of  accuracy 
may  be  encountered  in  shore  areas  where  the  available  empirical  data  are  marginal  in  quantity 
and/or  quality.  This  applies  particularly  to  the  Arctic  and  Antarctic  coastlines,  where  large  ice 
sheets  may  also  cause  some  tidal  distortions  and  retardations  which  have  been  assumed  as 
negligible. 

(0  Interactions  of  Oceanic  and  Terrestrial  Tides:  All  N'SWC  partial  tide  models  include 
effects  of  tide-generated  terrestrial  and  oceanic  mass  perturbations  in  the  simplified  forms  des¬ 
cribed  in  Section  1.  In  agreement  with  other  investigations,  the  inclusion  of  the  earth  tides 
produced  significant  first-order  effects  on  the  tidal  amplitudes.  As  expected  useful  second-order 
corrections  of  amplitudes  and  phases  were  registered  when  oceanic  load  effects  were  introduced, 
which  is  also  in  agreement  with  other  numerical  experiments  conducted  for  instance  by  Zahel 
(1477).  Accad  and  Pekeris  (1478).  L-stes  (1479).  and  Parke  and  Hendershott  ( 1979). 

(g)  Hydrodynamical  Properties:  The  ocean  tidal  equations  are  the  basic  LTlis  augmented 
by  physically  meaningful  lateral  eddy  dissipation  and  linear  bottom  friction  terms  with  novel 
eddy  and  friction  coefficients  depending  linearly  on  the  respective  vertical  and  horizontal  cell 
wall  areas.  Accordingly,  eddy  dissipation  was  found  significant  in  deep  oceans  while  bottom 
friction  dominated  in  shallow  waters.  It  may  be  emphasized  that  both  terms  were  needed  to 
achieve  realistic  decay  and  dispersion  features  and  to  bring  the  computed  data  in  good  agreement 
with  empirical  data.  They  were  not  necessary  to  enforce  stability  of  the  numerical  procedure.  In 
particular,  the  friction  coefficients  were  meshwise  modified  by  the  unique  hydrodynamical  inter¬ 
polation  in  order  to  achieve  a  realistically  smooth  acceptance  of  over  2000  empirical  tide  data 
pairs. 

(It)  Mathematical  Features:  The  finite  difference  analogue  of  the  ocean  tidal  equations 
uses  the  Hansen-Zahel  (Zahel  1470)  staggered  difference  scheme  in  space  and  improved  differences 
in  time.  At  boundaries,  the  free-slip  condition  is  strictly  imposed  while  the  no-crossflow  condition 
is  acceptably  violated  by  the  hydrodynamical  interpolation  of  empirical  boundary  data,  l  o  avoid 
the  polar  singularity  of  the  tidal  equations,  a  new  series  expansion  is  incorporated,  which  will  be 
described  in  a  forthcoming  paper  (Schwiderski  1 98 1  i). 


V  APPLICATIONS  Ol  THL  NSWC  OCEAN  TIDE  DATA 


Various  applications  of  the  NSWC  ocean  tide  models  require  indispensable  consideration  ol 
the  major  features  summarized  in  Section  2.  This  is  particularly  true  in  applications  to  compute 
interactions  of  oceanic  and  terrestrial  tides,  which  manifest  themselves  as  anomalies  in  the  earth's 
gravity  field.  Though  such  perturbations  amount  to  only  a  few  microgals.  they  can  now  be  de¬ 
termined  by  precision  measurements  such  as  those  carried  out  by  Kuo  et  al.  (1970)  and  more 
recently  on  a  global  scale  by  Melchior  et  al.  I  1980).  On  the  other  hand,  gravity  anomalies  caused 
by  ocean  tidal  loading  can  be  computed  from  accurate  ocean  tide  (e.g..  Schwiderski  1980b  and 
108 lj)  and  solid  earth  te.g..  Tarrell  l')72)  models  (see,  for  instance.  (load  1080a.  band  Melchior 
et  al.  1080).  Unfortunately,  both  the  experimentally  derived  and  computed  gravity  perturbations 
of  ocean  tidal  loads  are  subject  to  the  following  difficulties. 

(a)  Empirical  Gravity  Anomalies  of  Ocean  Tidal  Loads.  The  measured  gravity  anomalies 
include  dominant  perturbations  caused  by  the  direct  terrestrial  tides,  which  are  usually  modeled 
by  neglecting  all  lateral  density  and  elasticity  variations  of  the  solid  earth.  Yet,  those  heteroge¬ 
neities  exist  and  may  not  be  negligible  especially  in  tectonically  active  areas.  Furthermore,  after 
subtraction  of  the  modeled  terrestrial  tide  effects,  the  remaining  gravity  anomalies  must  be  sub¬ 
jected  to  the  same  harmonic  analysis  as  the  measured  ocean  tide,  which  is  particularly  controversial 
for  short-time  observations. 

(b)  Computed  Gravity  Anomalies  of  Ocean  Tidal  Loads:  The  computation  of  the  yielding 
of  the  solid  earth  under  the  oceanic  tide  is  also  based  on  a  highly  simplified  lateral  homogeneous 
earth  model  as  the  terrestrial  tide  model  (a)  above.  Moreover,  as  has  been  pointed  out  by’various 
investigators  (e.g..  Pekeris  1978.  Acead  and  Pekeris  1978.  Goad  1980a.  b,  and  Schwiderski  1980a). 
the  purely  numerical  computation  of  gravity  anomalies  from  ocean  tidal  loads  poses  considerable 
difficulties  and  may  be  the  cause  of  significant  errors.  This  is  particularly  true  at  coastal  and 
oceanic  stations  where  the  convolution  of  the  numerical  ocean  tide  with  a  singular  Green's 
function  must  be  computed.  Indeed,  depending  on  the  actual  numerical  procedure  chosen,  the 
singularity  may  amplify  in  various  ways  any  or  ail  of  the  deficiencies  of  the  NSWC  ocean  tide 
models  discussed  in  Section  2. 

In  spite  of  all  the  possible  error  sources.  Goad  (1980a.  b)  and  Melchior  et  al.  (1980)  arrive 
at  excellent  agreements  for  the  bulk  of  their  compared  empirical  and  computed  gravity  anomalies. 
These  results  lend  convincing  support  to  the  accuracy  claim  of  Melchior  et  al.  concerning  their 
worldwide  gravity  measurements.  Vice  versa,  these  investigations  verify  indirectly  the  oceanwide 
accuracy  and  usefulness  of  the  NSWC  tide  models.  In  fact,  a  detailed  analysis  of  the  massive 
Melchior  et  al.  data  indicates  that  all  of  the  remaining  discrepancies  between  the  empirical  and 
computed  data  can  probably  be  traced  to  one  or  more  of  the  modeling  and  computing  difficulties 
elaborated  above. 


b 


<c)  Immediate  Improvements:  In  coastal  and  oceanic  areas  some  of  the  ilel'ieiencies  ol  the 
NSWC  tide  models  (Section  2,  h.  c.  e)  can  he  immediately  removed  In  preparing  more  detailed 
limited-area  tide  models  on  a  locally  refined  grid  system.  Models  of  this  sort  have  alreadx  hern 
constructed  by  various  researchers  especially  for  excluded  or  insufficiently  resolved  bordering 
waters.  Where  no  refined  models  are  available,  ihe  refinement  may  be  accomplished  by  linear  ot 
higher-order  interpolation  of  the  listed  1°  x  1°  data  and  by  using  additional  empirical  data.  say. 
from  the  huge  collections  of  the  British  Admiralty  ( 1  c> 7 7 )  and  the  International  Hydrographic 
Bureau  (1978),  which  are  in  possession  of  the  author.  In  such  a  refining  process  the  hydrody 
namieally  interpolated  empirical  tide  data,  which  have  been  displaced  to  the  centers  of  the  grid 
cells  (Section  2.  c>.  should  be  first  returned  as  closely  as  possible  to  the  original  tide  gauge 
locations.  Specifically,  in  boundary  cells  the  empirial  data  should  be  moved  from  the  cell  centers 
to  the  centers  of  the  boundary  segments  (f  igure  la).  For  small  oceanic  (unresolved)  islands,  the 
tabulated  empirical  data  should  be  moved  from  the  cell  centers  to  the  original  stations  (Figure  lb). 
This  relatively  simple  local  procedure  is  strongly  recommended  where  the  ocean  tides  sutler 
drastic  distortions  and  retardations  from  the  deep  oceans  to  the  shallow  areas. 


/ 

/ 

'I 

A'  * 

/ 

a.  In  Boundary  Cells 


b.  Near  Unresolved  Islands 


Figure  1 .  Replacing  of  Empirical  Tide  Data 


The  important  effect  of  the  replacement  of  empirical  data  may  be  illustrated  by  the  incredibly 
sharp  decay  of  the  tidal  ranges  (=  double  amplitudes)  across  the  continental  shelf  of  the  New  York 
Bight.  As  can  be  seen  in  Figure  2,  the  tidal  ranges  measured  by  Gill  and  Porter  (1980)  drop  from 
4.6  ft  at  the  coast  over  a  shelf  width  of  less  than  200  km  to  an  almost  constant  height  of  about 
3  ft  in  the  deep  ocean.  The  global  NSWC  ocean  tide  model  approximates  this  distortion  very 
closely,  but  only  when  the  empirical  boundary  datum  "x  is  returned  to  the  original  boundary 
position  x.  Obviously,  simple  linear  interpolations  between  the  adjusted  model  data  x  and  x' 
could  significantly  improve  computations  of  corresponding  gravity  anomalies  in  that  neighbor¬ 
hood. 
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Figure  2.  Sharp  Decay  of  Tidal  Ranges  Across  the  Continental 
Shelf  of  New  York  Bight.  ©Gill-Porter  Measurements,  (x,  x'l 
NSWC  Tide  Model,  x  Replaced  Empirical  Boundary  Datum  x) 


(d)  Future  Improvements:  Refined  local  tide  models  are  presently  under  consideration  by 
the  author.  The  results  of  those  models  will  subsequently  be  used  in  connection  with  forthcoming 
new  and  better  empirical  data  to  improve  the  global  NSWC  tide  models.  An  error  reduction  from 
10  cm  to  3  cm  is  already  being  called  for  in  various  applications.  Of  course,  a  model  of  this 
accuracy  will  require  more  and  better  measured  and  analyzed  empirical  data  especially  in  critical 
areas. 
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4.  NSWC  GLOBAL  OCLAN  TIDE  DATA  (GOTD)  TAPE 


The  NSWC  GOTD-1981  tape  contains  the  amplitudes  £'m  (in  meters)  and  Greenwich  phases 
51  (in  degrees,  0°  =  360°)  of  the  leading  harmonic  ocean  tidal  modes  (M2 .  S2 .  K, ,  (), .  N, .  I', 
K2.Q,,  Mf.  Mm,  Ssa)  on  a  1°  x  1°  spherical  grid  system.  This  GOTD-1981  tape  is  identical  to  the 
earlier  GOT1)  - 1 080  tape,  which  did  not  contain  the  Mm  and  Ssa  data. 

The  entire  data  set  (£‘  .  )  is  arranged  by  modes  according  to  the  mode  nttmbei 

i  =  1,  2.  3 .  11  as  defined  in  Table  2.  which  also  lists  the  corresponding  species  number  v  (com¬ 

pare  also  Table  1 ). 


Table  2.  Entire  Data  Set  (see  Table  1 ) 


Tidal  Symbol 


Mode  Number:  i 


M, 


S2  K, 
2  3 


O, 


N, 


K2  0,  Ml'  Mm  Ssa 
7  8b  10  11 


Species  Number:  v 


1 


10  0  1) 


In  each  mode  (i  =  1,2 . 11),  the  data  are  arranged  by  colatitude  numbers  n  =  1.2 . lox 

(see  Equation  2)  in  consecutive  pairs  of  blocks  each  of  which  contain  361  words  of  amplitudes 

or  phases  arranged  by  longitude  numbers  m  =  1 , 2 .  360  (see  Equation  1 )  with  the  last  (361  -st ) 

word  indicating  the  colatitude  number  n.  The  first  block  in  each  pair  contains  the  amplitudes 
^  n  in  Format  F5.3  and  the  second  one  the  phases  Sjn  n  in  Format  F5.1.  The  data  have  been 
bloekwise  generated  by  the  BUFFER-OUT  statement  on  the  CDC  6700  computer.  The  magnetic 
tape  has  the  following  properties:  seven  track,  BCD  form,  even  parity,  556  bpi.  and  unlabeled. 

The  tidal  constants  (£jn  5jn  n)  are  representative  for  the  1°  x  1°  square  with  the  geo¬ 
graphical  center  point  (X|n  ,  0  ).  where 

Xm  =  (m  -  0.5)°  =  longitude  east  (m  =  1.2 . 360)  ( I  » 

0 n  =  (n  -  0.5)°  =  colatitude  (n  =  1,2 . 168)  (2) 

On  land  (continental  or  large  islands;  see  Scliw iderski  1978c)  all  tidal  constants  have  been  set  to; 

£in  n  =  9-999-  „  =  999.9  for  land  (3 1 


9 


The  precise  definition  of  those  data,  their  accuracy  and  limitation,  and  their  proper  applica¬ 
tion  and  improvement  lias  been  discussed  in  Sections  2  and  3.  To  avoid  possible  misleading  con¬ 
clusions  any  application  of  those  data  must  evaluate  realistically  the  listed  features  and  suggestions. 


The  instantaneous  ocean  tide  (in  meters)  of  any  model  i  (=  1,  2, ,  1 1)  can  be  found  at 
any  given  instant  (y  =  year,  d  =  day  number  of  year  y.  and  t  =  universal  time  of  day  d  in  seconds) 
and  any  given  ocean  point  (X,  0)  by  linearly  interpolating  the  amplitude  and  phase  6j  belonging 
to  the  point  (X,  0)  in  the  corresponding  partial  tide  tables  and  computing 

(X,  0.  t)  =  cos  [o|  t  +  7r  (Xj  -  Sj)/ 1803  (4) 

In  this  formula  <r  and  xt  denote,  respectively,  the  frequency  and  the  astronomical  argument  of 
the  partial  tide  listed  in  Table  1 . 


The  corresponding  instantaneous  astronomical  partial  equilibrium  tide  ^  (or  tide-generating 
potential  CT;  see  Schwiderski  1978a,  1979c,  d,  or  1980a)  is  determined  by 


=  b,  (X,d.  t)  =  f,(  (0)  cos  [o.  t  +7r  (Xj  +  v.  X)/180] 


(5) 


where  K.  =  K.  and  v  =  v  denote,  respectively,  the  tidal  amplitude  and  species  number  listed  in 
Table  1 .  The  species  function  of  colatitude  0  is  defined  by 


(„{0)=  | 


sin2  0  for  v  -  2 
sin2  0  for  v  -  1 

—  (3  sin2  0-2)  for  v  -  0 


(6) 


In  a  simple  second-order  Love-number  approximation,  the  corresponding  solid-earth  tide 
is  given  by  (see  Schwiderski  1978a,  1979c,  d,  or  1980a) 

r;  =  r;  (A.  (?.t)^o.6i  2^  <?> 

Similarly,  in  a  simple  Accad-Pekeris  approximation  (see  Schwiderski  1978a.  1979c.  d,  or  1980a), 
the  corresponding  dip  (yielding)  of  the  solid  earth  in  response  to  the  oceanic  tidal  load  f  is 
given  by 

(X.  0.  t)  ^  -  0.0667  (8) 

Thus,  the  corresponding  instantaneous  geocentric  partial  tide  (relative  to  the  geoidal  sea 
surface  at  rest)  is  determined  by  linear  superposition  to 

r?  =  tv(x.<u)  =  r.  +  £  +  {f°  (9) 
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Finally,  by  linear  superposition  of  all  corresponding  partial  tides,  one  finds  the  instantaneous 
local  tides: 


f  =  f  (X,  0,  t)  =  =  total  oceanic  tide  (10a) 

i 

=  f1  (X,  0,  t)  =  E  ff  =  total  geocentric  tide  (10b) 

i 

Note  to  Users:  Copies  of  the  NSWC  GOTD  tape  in  its  discribed  standard  format  (or  any  other 
desired  format)  may  be  requested  from  the  authors: 


F.  W.  Schwiderski,  Code  K104,  Tel.  (703)  663-8406 
L.  T.  Szeto,  Code  K 14,  Tel.  (7 03)  663-7741 
NSWC,  Dahlgren,  Virginia  22448 
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5.  NSWC  RANDOM-POINT  TIDE  PROGRAM  (RPT1DE) 


The  following  NSWC  RPT1DE  computes  from  the  harmonic  tidal  constants  listed  on  the 
NSWC  GOTD  1981  (or  GOTD  1980)  tape  (Section  4)  the  instaneous  oceanic  and/or  geocentric 
tides  at  given  random  points  and  instances.  The  corresponding  Users’  Guide  and  Program  Listing 
are  given  in  Appendixes  A  and  B.  A  fast  program  to  compute  geocentric  tides  along  satellite  tracks 
is  described  in  Sehwiderski,  E.  W.  and  Szeto,  L.  T.  1981. 

1)  Input  Data 

(1)  (£'  ,5'  n)  =  ocean  tide  amplitudes  (in  m)  and  Greenwich  phases  (in  deg)  from 

GOTD  1981  (or  1980)  tape,  where 

i  =  1 , 2,  3, ...  ,  1 1  =  mode  numbers  for  GOTD  1981 
(0  <  i  <  10  for  GOTD  1980) 

m  =  1,2, ...  ,  360  (=0)  =  longitude  numbers 

n=  1,2, ...  ,  168  =  colatitude  numbers 

Note:  £'  =  9.999,  5‘  n  =  999.9  on  land 

,n  *  m  ,n 

6‘m  n  =  360°  =  0°  (phase  jump) 

(2)  y  =  year  >  1975  (fixed  for  one  run!) 

(3)  d  =  day  of  year  y  (d  =  1  for  January  1st,  also  fixed  for  one  run!) 

(4)  tj  =  specified  instances  (in  sec  after  Greenwich  midnight  of  day  d,  where  j  =  1,  2,  3, ...  ,J 

(5)  (Xj,  0.)  =  longitudes  (East)  and  latitudes  (+  North,  -  South)  in  deg.  of  random  points 

corresponding  to  instances  t. 

(6)  J  =  number  of  random  points 

(7)  I  =  number  of  superposed  tidal  modes  if  1  <  l  <  1 1  with  option 

-  1  <  I  <  -  1 1  for  computation  of  single  tidal  mode  i  =  |1|,  |1|  <  10 
for  GOTD  1 980  tape) 

(8)  a  =  0  or  =  1  =  earth  tide  parameter  ( a  =  1  earth  tide  included,  a  =  0  earth  tide  excluded) 
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(9)  (3  =  0  or  =  1  =  ocean  load  parameter  </3  =  1  ocean  load  included,  0  =  0  ocean  load  excluded) 


(10)  7  =  0  or  =  1  =  ocean  tide  parameter  (7  =  1  ocean  tide  included,  7  =  0  ocean  tide  excluded) 


II)  Preliminary  Computations 


With  the  tidal  parameters  v  -  v.,  a  -  a.,  K  =  Kj;  ho,  so,  po  specified  in  Tables  1  and  2  compute: 


h  =  ir  hQ  / 1 80, 


and  set  (Et  =  0.612  1C) 

E,  =0.148  308. 
E4  =0.061  515, 
E7  =  0.018  791, 
E,  o  =0.013  480, 


s  =  ^so/180,  p  =  7rpo/180 
X  =  X,  for  i  =  1,  2,  3, ....  1 1 


E2  =0.069  059, 
Es  =  0.028  396, 
E8  =  0.01 1  785, 
E,  =0.011  901, 


E3  =0.086  638, 
E6  =  0.028  668, 
E9  =  0.025  546, 


also 


7  =  7-  0.0667  0 


III)  Main  Computation 

(1)  Set  j  =  1 

(2)  Set  fj  =  0  and 

=  1  for  I  >  0, 
i 

=  |I|  for  1  >  0  (single  mode!) 


(3)  Compute: 


m  =  Int  [Xj  +  0.5  J  +  1 ,  0  =  m  -  (X.  +  0.5) 
n  =  Int  (90.5  -0j|  +  l,0  =  fl  =  n-  (90.5  -  0.) 


and  transfer  to  core  memory 


-  I  ,  n-  1  ’  -  1 ,  n -  1  ,  n-  1  ’  ,  n-  1  ^ 


(pi  8‘  )• 

I  ,  n  ’  um  -  I  .  n  ’’ 


(*) 
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where  m  -  1  =  0  -*•  360 


Check  for  land  points  and  replace 

(a)  if  tJn  n_  ]  =  9.999,  replace  0  -*  0 

(b)  if  ^  n  =  9,999,  replace  0  -*  1 

(c)  if  (a)  and  (b)  hold,  replace  \p  -*  1 

(d)  if  £'m  _  j  n_  ,  =  9.999.  replace  0  -*  0 

(e)  if  {‘m  _  j  n  =  9.999,  replace  0  1 

(f)  if  (d)  and  (e)  hold,  replace  \p  -*  0 

(4)  Interpolate  or  the  ruled  second-order  surface  in  4>  and  6 : 

*j  =  ( 

+  n,i  -  0>€*m_i.  nl 


Test  for  360°  phase  jumps  and  replace 


if  8L 


if  8 


\ 

m  -  1 


n  -  1 


if  5*  .  -  6* 

m  -  1  ,  n  m  ,  n 


if8‘ 


m  -  I  .  n  -  1 


m  .  n  -  1 


^  *m- I  ,  n  ~  ®m.  n  -  I 


>  180,  replace  5'm  n_,  -  8*m  -  360 

<  -  1 80,  replace  5jn  n  -^5>m  n-  360 

>  1 80,  replace  8'ni_l  n_,  -360 

<  -180,  replace  5^  n  -►  S'm  n  -  360 

>  1 80,  replace  8*m  _  t  n  ~  360 

<  -  180,  replace  S1  -*■  5*  -  360 

>  180,  replace  5'  .  .  ~*8'  .  ,  -  360 

<  -180,  replace  6'  „  .  ->  6‘  .  -  360 

1  m  ,  n -  I  m  *  n~  I 

>180.replace6|n_lin-fiini_1>n-360 

<  "  180,  replace  8;ii  n_,  -  S'm  -  360 


(**) 
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if  S'  ,  -  5* 

m  -  1  ,  n  in  -  1  .  n  -  1 


>  180.  replace  5;n_ ^  n  ->6;n_K  n 


360 


<  -180,  replace  -  S', .  ,  ,  -  360 


Use  adjusted  5s  to  interpolate  S'  by  formula  (**)  with  £  replaced  by  6. 


Now,  compute  and  replace 


fj  -*■  +  7  £|  eos  [a;  f  +  X;  -  ir  6‘/1801 
+  a  (<t>j )  cos  [oi  tj  +  x,  +  7r  v.  Xj/180] 


where 


E  cos2  0  for  v  =  2 

i  j  i 

E  sin  2  0,  for  v.  =  1 

i  j  i 

l 


v  -  Ej  (3  cos2  0^  -  2)  for  =  0 


If  i  <  |I|,  replace  i  -*  i  +  1 .  transfer  to  core  memory  new  data  (*),  and  repeat  routine  (4). 

If  i  <  |1|,  print  A.,  0.,  tj .  and  fj- 

If  j  <  J,  replace  j  -*j  +  1  and  repeat  (2)  through  (4). 

If  j  <  J  stop  program. 


Note:  This  program  does  not  consider  any  special  features  or  suggestions  made  in  Sections  2 
and  3.  Hence,  if  a  higher  accuracy  near  coastal  boundaries  is  needed,  an  appropriate  local  compu¬ 
tation  should  be  used  following  the  suggestions  of  Section  2. 
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APPENDIX  A 


RPT1DE  PROGRAM  USER’S  GUIDE 

Fho  Random-Point  Tide  <  RP 1 11)1  )  program  computes  instantaneous  oceanic  and/or 
geocentric  tides  at  randomly  given  geographical  points  and  specified  instances  as  described  in 
Section  5.  The  program  listing  in  Appendix  B  shows  RPTIDL  written  in  FORTRAN  IV  extended 
for  the  CDC  6700  computer  under  the  SCOPE.  3.4  operating  system.  The  BUFFER  IN  statement 
buffers  in  GOTD  1980  or  1981  tide  constants  (Section  4).  which  are  in  binary  mode.  The  GOTO 
magnetic  tape  is  in  coiled  mode:  its  files  are  attached  as  follows  for  the  example  I  =  6.  L  =  I  I  I 
(see  Sections  4  and  5 ): 


attach,  tape:  l  z'tn  n 

ATTACH.  TAPE  L+  1.  hlm  n 
ATTACH.  TAPE  L  +  2.  n 
ATTACH.  TAPE  L  +  3.  6^  n 

ATTACH.  TAPE  L  +10.  n 

ATTACH.  TAPE  L+  11.  8*  „ 

All  other  input  data  are  entered  using  data  cards  consisting  of  one  type  A  data  card  one  or  more 
type  B  data  cards.  The  type  A  card  contains  8  right-justified  integers  with  the  format  of  815. 
The  integers  represent  the  following  ordered  variables  (Section  5  symbols  in  parentheses): 

YEAR  =  (y )  =  year  >  1975  (fixed  for  one  run!) 

DAY  =  (d)  =  day  of  YEAR  (d  =  1  for  January  1st.  also  fixed  for  one  run!) 

JPTS  =  ( J )  =  total  number  of  random  points  or  instances  with  same  YEAR  and  DAY 

MTIDE  =  (1)  =  total  number  of  tidal  modes  to  be  superposed  (1  <  i  <  1  where  I  <  10  for  GOTD 
1980  and  1  <  12  for  GOTD  1981).  Option  if  MTIDE  =  -  N  with  0  <  N  <  10  or 
1 2,  the  computed  tide  is  the  single  partial  tide  i  =  N 

ETINC  =  (a)  =  earth  tide  parameter  (ET1NC  =  0  excludes  and  ETINC  =  1  includes  earth  tide) 

LOAD  =  (0)  =  ocean  load  parameter  (LOAD  =  0  excludes  and  LOAD  =  1  includes  ocean  load) 


/*■ 


- 1 .  ir  Hi  I  II ■  ■uni 


OCEAN  =  (7)  =  ocean  title  parameter  (OCEAN  =  0  excludes  anti  OCEAN  =  I  includes  ocean  tide) 


IT  l  =  beginning  unit  number  of  2  *  :MT1DE|  consecutive  unit  numbers  to  which  GOT!)  tiles 
are  attached  (see  example  above) 

Each  type  B  data  card  contains  three  numbers  with  the  format  3F20.8.  These  numbers  are, 
respectively: 

TIMM  =  (t)  =  time  (in  see)  of  given  instant  after  Greenwich  midnight  of  day  DAY  (may  be  several 
days  later) 

LONG  =  (X)  =  Longitude  (Last:  in  degrees) 

LAT  =  (<p)  =  Latitude  (+  North,  -  South;  in  degrees) 

Other  major  notations  corresponding  to  Section  5  are: 


I)  =  (6)  =  Oceanic  tidal  phase 
X  =  (£)  =  Oceanic  tidal  amplitude 
T  =  (f)  =  instantaneous  tidal  height 
NU  =  (n)  =  tidal  species  number 
FRLQ  =  (ct)  =  tidal  frequeaney 
ASTRO  =  (x)  =  astronomical  argument 
LT1DL  =  (L)  =  earth  tide  constant  amplitude 
LT  =  (H(0))  =  earth  tide  geographical  amplitude 

Note:  For  land  points  in  squares  with  9.999  vertex  amplitudes  only  the  earth  tide  is  meaningful 
if  ETINC  =  1 ,  LOAD  =  OCEAN  =  0.  For  improved  accuracy  in  coastal  waters,  special 
computations  should  incorporate  the  suggestions  of  Section  3  c  and  d. 
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RPTIDE  PROGRAM  LISTING 
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SUBROUTINE  CO"f>UT 
(  M  t  X,  0,  XN»  ON»  LONG 
,  FUi  TINE,  ASTRO.  S.L.ETTNC  .GAHHAR 
.  OTT.  OU.  ETT 
.  FPFO.  LAT,  ETIOF»  Ti 


OEFINEO  1 

IAT  REH  F.P.  RFF  S  9  17  55  2*58  DEFINED  1 

LONG  PE»L  F.P.  RE  F  5  9  16  65  DEFINED  1 

P  INTEGER  F.P.  REFS  12  13  16  19  20  2*23  2*27 

29  30  2*33  35  35  36  37 
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